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Abstract 

The confinement mechanism proposed earher by the author is employed for to com- 
pute the decay constants fp corresponding to leptonic decays P ^ l^ + 1^1, I = fj,,e, 
where P stands for any meson from vr^, K^. For this aim the weak axial form 
factor of P-meson is nonperturbatively calculated. The study entails estimates for 
parameters of the confining SU(3)-gluonic field in charged pions and kaons. The 
corresponding estimates of the gluon concentrations, electric and magnetic colour 
field strengths are also adduced for the mentioned field at the scales of the mesons 
under consideration. Further the obtained results are applied to the problem of chi- 
ral symmetry breaking in quantum chromodynamics (QCD). It is shown that in 
chirally symmetric world masses of pions and kaons are fully determined by the 
confining SU(3)-gluonic field among (massless) u, d and s quarks and not equal 
to zero. Accordingly chiral symmetry is sufficiently rough approximate one holding 
true only when neglecting the mentioned SU(3)-gluonic field between quarks and no 
additional mechanism of the spontaneous chiral symmetry breaking connected to 
the so-called Goldstone bosons is required. Finally, a possible relation of the results 
obtained with a phenomenological string-like picture of confinement is discussed 
too. 
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1 Introduction 



In Refs. [1-3] for the Dirac- Yang-Mills system derived from QCD-Lagrangian 
an unique family of compatible nonperturbative solutions was found and ex- 
plored, which could pretend to describing confinement of two quarks. The 



Preprint submitted to Elsevier 



10 September 2008 



applications of the family to the description of both the heavy quarkonia 
spectra [4,5] and a number of properties of pions, kaons, t]- and ^^'-mesons [6- 
8] showed that the confinement mechanism is qualitatively the same for both 
light mesons and heavy quarkonia. At this moment it can be described in the 
following way. 

Two main physical reasons underlie linear confinement in the mechanism im- 
der discussion. The first one is that gluon exchange between quarks is realized 
with the propagator different from the photon-like one and existence and form 
of such a propagator is a direct consequence of the unique confining nonpertur- 
bative solutions of the Yang-Mills equations [2,3]. The second reason is that, 
owing to the structure of mentioned propagator, quarks mainly emit and inter- 
change the soft gluons so the gluon condensate (a classical gluon field) between 
quarks basically consists of soft gluons (for more details see Refs. [2,3]) but, 
because of that any gluon also emits gluons (still softer), the corresponding 
gluon concentrations rapidly become huge and form a linear confining mag- 
netic colour field of enormous strengths, which leads to the confinement of 
quarks. This is by virtue of the fact that just magnetic part of the mentioned 
propagator is responsible for a larger portion of gluon concentrations at large 
distances since the magnetic part has stronger infrared singularities than the 
electric one. Under the circumstances physically nonlinearity of the Yang-Mills 
equations effectively vanishes so the latter possess the unique nonperturbative 
confining sohitions of the Abelian-like form (with the values in Cartan subal- 
gebra of SU(3)-Lie algebra) [2,3] which describe the gluon condensate under 
consideration. Moreover, since the overwhelming majority of gluons is soft 
they cannot leave hadron (meson) until some gluons obtain additional energy 
(due to an external reason) to rush out. So we also deal with confinement of 
gluons. 

The approach under discussion equips us with the explicit wave functions for 
every two quarks (meson or quarkonium). The wave functions are parametrized 
by a set of real constants aj,bj, Bj describing the mentioned nonperturbative 
confining SU(3)-gluonic field (the gluon condensate), and they are nonpertur- 
bative modulo square integrable solutions of the Dirac equation in the above 
confining SU(3)-field and also depend on fj,o, the reduced mass of the current 
masses of quarks forming meson. It is clear that under the given approach just 
constants aj,bj, Bj, fiQ determine all properties of any meson (quarkonium), i. 
e., the approach directly appeals to quark and gluonic degrees of freedom as 
should be according to the first principles of QCD. Also it is clear that the 
constants mentioned should be extracted from experimental data. 

Such a program has been to a certain extent advanced in Refs. [4-8]. The 
aim of the present paper is to continue obtaining estimates for aj,bj,Bj for 
concrete mesons starting from experimental data on spectroscopy of one or 
another meson. We here consider charged pions and kaons and their leptonic 
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decays P ^ l"^ + ui, I — ii,e, where P stands for any meson from tt^, K"^. 

Of course, when conducting our considerations we shaU rely on the standard 
quark model (SQM) based on SU(3)-flavor symmetry (see, e. g., Ref. [9]), so 
in accordance with SQM tt^ = ud, ud, — us, us respectively. 

Under the circumstances Section 2 contains main relations underlying our 
approach. Section 3 is devoted to computing the electric form factor, the root- 
mean-square radius < r > and the magnetic moment of the mesons under 
consideration in an explicit analytic form while Section 4 gives an explicit 
analytic expression of the weak axial form factor for those mesons. Results of 
Sections 3 and 4 are used in Section 5 for obtaining estimates for parameters 
of the confining SU(3)-gluonic field for the mesons under discussion. Section 
6 employs the obtained parameters of SU(3)-gluonic field to get the corre- 
sponding estimates for such characteristics of the mentioned field as gluon 
concentrations, electric and magnetic colour field strengths at the scales of 
the mesons in question. In Section 7 results of previous sections are applied to 
the problem of chiral symmetry breaking in QCD. Section 8 explores a pos- 
sible connection of the results obtained with a phenomenological string-like 
picture of confinement while Section 9 is devoted to discussion and concluding 
remarks. 

At last. Appendices A and B contain the detailed description of main building 
blocks for meson wave functions in the approach under discussion, respectively: 
eigenspinors of the Euclidean Dirac operator on 2-sphere §^ and radial parts 
for the modulo square integrable solutions of Dirac equation in the confining 
SU(3)-Yang-Mills field. 

Further we shall deal with the metric of the fiat Minkowski spacetime M that 
we write down (using the ordinary set of local spherical coordinates r, i), ip for 
the spatial part) in the form 

ds^ = g^^dx^" dx'' = de - dr^ - r^d^"^ + sin^ Mip^) . (1) 

Besides, we have \6\ = \ det(5f^j,)| = (r^sint?)^ and 0<r<oo, 0<'i?<7r, 
< < 27r. 

Throughout the paper we employ the Hcaviside-Lorentz system of units with 
h = c = 1, unless explicitly stated otherwise, so the gauge coupling constant g 
and the strong coupling constant ag are connected by relation g^/i^in) = ag. 
Further, we shall denote by L2{F) the set of the modulo square integrable com- 
plex functions on any manifold F furnished with an integration measure, then 
L2{F) will be the n-fold direct product of L2{F) endowed with the obvious 
scalar product while f and * stand, respectively, for Hcrmitian and complex 
conjugation. Our choice of Dirac 7-matrices conforms to the so-called standard 
representation and is the same as in Ref. [6]. At last ® means the tensorial 
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product of matrices and /„ is the unit n x n matrix so that, e.g., we have 

\ 

/g = 7^ 

V 07^/ 

for any Dirac 7-matrix 7^* and so forth. 

When calculating we apply the relations 1 GcV~^ ~ 0.1973269679 fm, 1 ^ 
0.658211915 X 10-^4 GeV , 1 V/m ^ 0.2309956375 x lO'^^ GeV^ 1 T = 
Att X 10-^H/m X 1 A/m « 0.6925075988 x lO'^^ GeVl 

Finally, for the necessary estimates we shall employ the Too-component (volu- 
metric energy density ) of the energy-momentum tensor for a SU (3) -Yang-Mills 
field which should be written in the chosen system of units in the form 

T,. = + If;^ F:,g^^g'''g,^ . (2) 



2 Main relations 



As was mentioned above, our considerations shall be based on the unique fam- 
ily of compatible nonperturbative solutions for the Dirac- Yang-Mills system 
(derived from QCD-Lagrangian) studied at the whole length in Refs. [1-3]. Re- 
ferring for more details to those references, let us briefly describe and specify 
only the relations necessary to us in the present paper. 

One part of the mentioned family is presented by the unique nonperturbative 
confining solution of the SU (3)- Yang-Mills equations for gluonic field A — 
Af^dx^ = A^^Xadx^^ (A(j are the known Gell-Mann matrices, // = t,r, "d, </?, 
a = 1, 8) and looks as follows 



Alt = Al + -^Al = + ,^2t = -Al + -^Al = -^ + A2 



a/3 

Ai^ = Al + -^Al = bir + Bi , A2^ = -A^ + -^Aj = + B2 

2 



= = -(^1 + ^2)r - (Bi + B2) (3) 

with the real constants aj, Aj, bj, Bj parametrizing the family. The word unique 
should be understood in the strict mathematical sense. In fact in Ref. [2] the 
following theorem was proved: 
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The unique exact spherically symmetric (nonperturbative) solutions (i.e. de- 
pending only on r) of SU(3)- Yang-Mills equations in Minkowski spacetime 
consist of the family of (3). 

It should be noted that solution (3) was found early in Ref. [1] but its unique- 
ness was proved just in Ref. [2] (see also Ref. [3]). Besides, in Ref. [2] (see also 
Rcf. [6]) it was shown that the above unique confining solutions (3) satisfy 
the so-called Wilson confinement criterion [10]. Up to now nobody contested 
this result so if we want to describe interaction between quarks by spheri- 
cally symmetric SU(3)-fields then they can be only the ones from the above 
theorem. 

As has been repeatedly explained in Refs. [2-4,6], parameters Ai^2 of solution 
(3) are inessential for physics in question and we can consider Ai = A2 — 0. 
Obviously we have Z^j=i Ajt — Z^f=i Aj^ — which reflects the fact that for 
any matrix T from SU(3)-Lie algebra we have TrT = 0. Also, as has been 
repeatedly discussed by us earlier (see, e. g., Refs. [2,3]), from the above form it 
is clear that the solution (3) is a configuration describing the electric Coulomb- 
like colour field (components A^'^) and the magnetic colour field linear in r 
(components A^^) and we wrote down the solution (3) in the combinations 
that are just needed further to insert into the Dirac equation (4). 

For the sake of completeness one should note that the similar unique confining 
solutions exist for all semisimple and non-semisimple compact Lie groups, in 
particular, for SU(A^) with > 2 and \J{N) with N >1 [2,3]. Explicit form of 
solutions, e.g., for SU(A") with N = 2,4 can be found in Ref. [3] but it should 
be emphasized that components linear in r always represent the magnetic 
colour field in all the mentioned solutions. 

Another part of the compatible nonperturbative solutions for the SU(3)- Dirac- 
Yang-Mills system is given by the unique nonperturbative modulo square in- 
tegrablc solutions of the Dirac equation in the confining SU(3)-field of (3) 
^ = (^1,^2! ^3) with the four-dimensional Dirac spinors representing 
the jth colour component of the meson, so \I' may describe relative motion 
(relativistic bound states) of two quarks in mesons and the mentioned Dirac 
equation looks as follows 





( Hi<^i 

H2^2 I , (4) 
V i^3^3 



where Hamiltonian Hj is 



r \ 2 J r smv \ 2 2 
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-9l"(-M,t+l'^-,A,^ (5) 

V rsmv J 

with the gauge couphng constant g while /iq is a mass parameter and one 
should consider it to be the reduced mass which is equal, e. g., for quarkonia, 
to half the current mass of quarks forming a quarkonium. 

Then the unique nonperturbativc modulo square intcgrable solutions of (4) 
are (with Pauli matrix cxi, for more details sec Rcfs. [1,3]) 

with the 2D eigenspinor $j = ^ ^"'^ ^ of the Euclidean Dirac operator Vq on 

the unit sphere while the coordinate r stands for the distance between 
quarks. The explicit form of $j is discussed in Appendix A. We can call the 
quantity cuj relative energy of jth colour component of meson (while ipj is wave 
function of a stationary state for jth colour component) but we can see that 
if we want to interpret (4) as equation for eigenvalues of the relative motion 
energy, i. e., to rewrite it in the form Hip = ujip with ip = {ipi,ip2, "^s) then we 
should put LU = LJj for any j so that Hjipj = Wjipj = uripj. Under this situation, 
if a meson is composed of quarks gi_2 with different flavours then the energy 
spectrum of the meson will be given by e = m^^ + m^j + cu with the current 
quark masses m^^ ( rest energies) of the corresponding quarks. On the other 
hand for determination of cuj the following quadratic equation can be obtained 
[1-3] 

[g^a'^j+{nj+ajf]uj^-2{Xj-gBj)g\bj ujj+[{Xj-gBj)'^-{nj+aj)'^]g%'^j-nl{nj+ajf = 0, 



(7) 



that yields (at 5* 7^ 0) 



Ajg\bj ± {rij + aj)^J{nj + 2njaj + Aj)/xg + g%j{n^j + 2njaj) . 

(8) 

where 03 = -(ai + as), h = -(^'i + h), B3 = -{Bi + B2), Aj = - gBj, 



aj = y A| — g'^aj, rij = 0, 1, 2, while Xj = ±(/j + 1) are the eigenvalues of 
Euclidean Dirac operator Vq on unit sphere with Ij = 0, 1,2, .... It should be 
noted that in the papers [1-4,6] we used the ansatz (6) with the factor e*'^^* 
instead of e""^^* but then the Dirac equation (4) would look as —idt^ — H^i 
and in equation (7) the second summand would have the plus sign while the 
first summand in numerator of (8) would have the minus sign. In the papers 
[5,7,8] we returned to the conventional form of writing Dirac equation and this 
slightly modified the equations (7)-(8). In the given paper we conform to the 
same prescription as in Refs. [5,7,8]. 
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In line with the above we should have u — ui = uj2 = ^3 in energy spectrum 
e = + rriq^ + uj for any meson (quarkonium) and this at once imposes two 
conditions on parameters aj, bj, Bj when choosing some experimental value for 
e at the given current quark masses m^^ , niq^ . 

The general form of the radial parts of (6) is considered in Appendix B. Within 
the given paper we need only of the radial parts of (6) at rij = (the ground 
state) that are [see (B.5)] 

with Pj = yJiJ,Q — LJj+ g^j while Cj is determined from the normalization con- 
dition J^{\Fji\'^+\Fj2\^)dr = |. Consequently, we shall gain that *j e Lj{R^) 
at any t & R and, as a result, the solutions of (6) may describe relativistic 
bound states (mesons) with the energy (mass) spectrum e. 



2.1 Nonrelativistic and the weak coupling limits 



It is useful to specify the nonrelativistic limit (when c — > oo) for spectrum 
(8). For that one should replace g — > g/VKc, aj — > aj/VKc, bj — > bjVkc, 



Bj/yhc and, expanding (8) in 2; = 1/c, we shall get 



±Aioc 



It 



g^a] 



+ 



ujj{nj,lj,Xj) 
Xjg^ajbj 



Hrij + \Xj\y 



9^Bja]f{nj, Aj) 



where /(n,-. A,) = 4:Xjnj{n] + X]) + ^ (n| + Gn^Af + A|) . 



z+0{z^) 
(10) 



As is seen from (10), at c — > 00 the contribution of linear magnetic colour 
field (parameters bj, Bj) to spectrum really vanishes and spectrum in essence 
becomes the purely nonrelativistic Coulomb one (modulo the rest energy). 
Also it is clear that when 



00, LOj 



At last, one should specify the weak coupling limit of (8), i.e., the case g' — > 0. 
As is not complicated to see from (8), cuj ±/xo when g 0. But then 

quantities Pj = J fi^ — cj^ + g'^b'j — > and wave functions of (9) cease to 



be the modulo square integrable ones at g = 0, i.e., they cease to describe 
relativistic bound states. Accordingly, this means that equation (8) does not 
make physical meaning at = 0. 
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We may seemingly use (8) with various combinations of signes (±) before sec- 
ond summand in numerators of (8) but, due to (10), it is reasonable to take 
all signs equal to plus which is our choice within the paper. Besides, as is not 
complicated to see, radial parts in nonrelativistic limit have the behaviour of 
form Fji,Fj2 ~ r'^"*"^, which allows one to call quantum number Ij angular 
momentum for j'th colour component though angular momentum is not con- 
served in the field (3) [1,3]. So for meson (quarkonium) under consideration 
we should put all Ij = 0- 



2.2 Eigenspinors with A = ±1 

Finally it should be noted that spectrum (8) is degenerated owing to degen- 
eracy of eigenvalues for the Euclidean Dirac operator "Dq on the unit sphere 
S^. Namely, each eigenvalue of Vq X = ±(/-|- 1),/ = 0,1, 2..., has multiplicity 
2(/ + 1) so we has 2(/ + 1) eigenspinors orthogonal to each other. Ad referen- 
dum we need eigenspinors corresponding to A = ±1 (Z = 0) so here is their 
exphcit form [see (A. 16)] 

A = -1 : ^ = ? f ) e^^/^ or $ = ^ f ] e'^^l^ 

2 ^6-^2 ) ' 2 V-e-*2 j 

A = 1 : $ = - f^? e^^/^ o.^ = -{ "C" \ e-^/^ (11) 
with the coefficient C — 1/\/2t^ (for more details see Appendix A). 

3 Electric form factor, the root-mean-square radius and magnetic 
moment 

Within the present paper we shall use relations (8) at rij = = Ij so energy 
(mass) of mesons under consideration is given by mp = rUq^ + rrig^ + uj with 
CO — u!j{0, 0, Xj) for any j = 1, 2, 3 whereas 

^ _ g'^aih _^ aijjLQ _ g'^a2h _^ 0:2/^0 _ g'^ash _^ QsA^o _ f rnT,± - ruu - ma, 
~ Ai |Ai| ~ A2 IA2I ~ A3 IA3I ~ \ mK±-mu-m,,,, 

(12) 

and, as a consequence, the corresponding meson (quarkonium) wave functions 
of (6) are represented by (9) and (11). 
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3. 1 Choice of quark masses and the gauge coupling constant 



It is evident for employing the above relations we have to assign some values 
to quark masses and gauge coupling constant g. In accordance with Ref. [9], 
at present the current quark masses necessary to us are restricted to intervals 
1.5 MeV < m„ < 3 MeV, 3.0 MeV < < 7 MeV, 95 MeV < < 120 MeV, 
so we take m„ = (1.5 + 3)/2 MeV = 2.25 MeV, = (3 + 7)/2 MeV = 
5 MeV, rris = (95 + 120)/2 MeV = 107.5 MeV. Under the circumstances, 
the reduced mass /Iq of (5) will be rriumdl {rriu + rud) or mums/{mu + mg). 
As to gauge coupling constant g = y^Anas, it should be noted that recently 
some attempts have been made to generalize the standard formula for ag = 
C(s{Q^) = 127r/[(33 — 2n/) In (Q^/A^)] (n/ is number of quark flavours) holding 
true at the momentum transfer — > oo to the whole interval < ^/Q^ < 
oo. We shall employ one such a generalization used in Refs. [11]. It looks as 
follows {x = VQ^ in GeV) 

/ X _ 127r fi{x) 

(33-2n/)ln5!±^ ^ ^ 



with 



-1 



/i(x) = l+|(2^^^J^g-^ln^-l) \o.6x'-' j , Mx) = m\l+2.8x' 

wherefrom one can conclude that q;^ — > tt = 3.1415... when x — > 0, i. e., 
^ ^ 27r = 6.2831.... We used (13) at m = 1 GeV, A = 0.234 GeV, Uf = 3, 
X = m^± = 139.56995 MeV or x = mK± = 493.677 MeV to obtain g = 
6.091309951 necessary for our further computations at the mass scale of tt^- 
mesons oi g = 5.301208569 at the mass scale of i^T^-mesons. 



3.2 Electric form factor 



For each meson (quarkonium) with the wave function = (^j) of (6) we can 
define electromagnetic current = "^(/s (g) -f'')'^ = (\E't\E', \E't(/3 ^ a)^) = 
(p, J), a = 7°7. Electric form factor f{K) is the Fourier transform of p 

f{K) = J ^^^e-'^^^'d'^x = E / ^j^iC-^^^ci^a; = Y.fj{K) = 

Yl / (I^Jil^ + \Fj2\'^)^j^j^-d^x, d^x = sin^drdMif (14) 
j=i J r 

with the momentum transfer K. At rij = = Ij, as is easily seen, for 

any spinor of (11) we have — l/(47r), so the integrand in (14) does 
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not depend on Lp and we can consider vector K to be directed along z- 
axis. Then Kr = Kr cosd and with the help of (9) and relations (see Ref. 

[12]): S^T'^-^e-^dT = r(a)p-«, Re a,p > 0, r^-^e-^ (rs(Jr)) ^ 



r(a)(K^+p2)-'^/2 f^^^l^^^^f ^1^/^11 Y Re « > -1, Rep > |ImX|, T{a + 
^ V cos (a arctan (A/p)) y ' ^ 

1) = ar(a), /(f e-^-^^'=°^''sini?di? = 2 sin {K r) / {Kr), we shall obtain 



j(^) ^ X^/iW = E sin [2a, arctan (7^/(2/5,))] 



wherefrom it is clear that f{K) is a function of K"^, as should be, and we can 
determine the root-mean-square radius of meson (quarkonium) in the form 



<r > -- 



V 2aj + 3a, + l 



When calculating (15) also the fact was used that by virtue of the normal- 
ization condition for wave functions we have C|[P,?(1 — ghj/fijY -\- (5^(1 + 
^6,//3,)2] = (2/3,)2«.+V[3r(2a, + l)]. 

On the other hand, we can directly calculate < r > in accordance with the 



standard quantum mechanics rules as < r >= J r^vl/t^^x = y / r'^'^^'^jd^x 

and the result will be the same as in (16). So we should not call < r > of 
(16) the charge radius of meson (quarkonium)- it is just the radius of meson 
(quarkonium) determined by the wave functions of (6) (at n, = = /,) with 
respect to strong interaction, i.e., radius of confinement. Now we should notice 
the expression (15) to depend on 3- vector K. To rewrite it in the form holding 
true for any 4- vector Q, let us remind that according to general considerations 
(see, e.g., Ref. [13]) the relation (15) should correspond to the so-called Breit 
frame where — [^jig^ fixing metric by (1)] so it is not complicated 

to rewrite (15) for arbitrary Q in the form 

3 3 /2« A2a,+i sin[2Q;,-arctan(,AQ2T/(2/3 ))] 
f(Q^) = V /.(g2) = V- H= — ^ (17) 

which passes on to (15) in the Breit frame. 
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3.3 Magnetic moment 



We can define the volumetric magnetic moment density by m = q{v x J)/2 = 
(AiyJz — zJy)i + {zJx — xJz)^ + {xJy — yJj;)k]/2 with the meson charge q and 
J = ^'^(/a ck)^. Using (6) we have in the exphcit form 

j=i ' i=i ' 

- E(^i^.2 - i^^i^.i)^^ (18) 

with Pauli matrices (Ti,2,3- Magnetic moment of meson (quarkonium) is M = 
fy m(Px, where V is volume of meson (quarkonium) (the ball of radius < r >). 
Then at nj = Ij = 0, as is seen from (9), (11), Fj\ = Fji,F*2 = -Fj2, 

^'j(J20'i^j = for any spinor of (11) which entails Jx = Jy = 0, i.e., rUz = 
while Jy mx,yd^x = because of turning the integral over ip to zero, which is 
easy to check. As a result, magnetic moments of mesons (quarkonia) with the 
wave functions of (6) (at Ij — 0) are equal to zero, as should be according to 
experimental data [9]. 

Though we can also evaluate magnetic form factor F[Q'^) of meson (quarko- 
nium) which is also a function of (see Refs. [7,8]) but the latter will not be 
used in the given paper so we shall not dwell upon it. 



4 Weak eixial form factor 

4 . 1 Preliminaries 

Let us now address to leptonic decays P ^ l"^ + ui, I = n,e. According to the 
standard theory of electroweak interaction the width of those decays is given 
by (see, e.g., Refs. [9,14]) 

with the Fermi coupling constant Gp — 1.16637 x 10~^ GeV~^, lepton mass mi, 
the charged P-meson mass mp, the corresponding clement V of the Cabibbo- 
Kobayashi-Maskawa mixing matrix {V = Vud = cos'dc for P = and 
V — Vus = sint^c for P = with the Cabibbo angle i)c ~ 13°), while 
the information about strong interaction of quarks in P-mesons is encoded 
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in a decay constant fp. Fnrthcr parametrization in the form fp = mp^ en- 
tails $ ^ 0.9364479961 at fp = U ^ 130.7 McV and $ ^ 0.3236934271 at 
fp = Jk ^ 159.8 MeV [9]. We can notice that the only invariant which $ 
might depend on is Q'^ — nip, i.e. we should find such a function TiCp') for 
that :r(g2 ^ ^2^) ~ 0.9354479961 or 0.3236934271 respectively. It is obvi- 
ous from physical point of view that T should be connected with electroweak 
properties of P-mesons. 

Let us now remark that, as is well known [9,14], the mentioned leptonic decays 
of charged pions and kaons are governed by the axial- vector part of the weak 
charged hadronic current. We can try to construct that part from the P- meson 
wave function of (6) and the wave function of vacuum state which may be 
chosen from the considerations that it should be similar to (6) in the form but 
it should not contain any explicit information about the parameters a^, 6j, Pj, 
describing the P-meson properties in our approach and connected directly 
to quark and gluonic degrees of freedom. This requirement can be satisfied by 
^0 of the form 

-40 («-*"'''•-'( .J = 1.2.3 (20) 

with the 2D eigenspinor = ( j of the Euchdean Dirac operator Dq on 

the unit sphere and some real constant > whose value will be fixed 

below (see Section 5) while the phase factor e"*'^^* is left only for convenience 
and has no influence on the subsequent calculations. 



^.2 Weak axial form factor 



To fix spinor $^ let us consider the axial-vector part of the weak charged 
hadronic current in the form = 'l'o(-^3 ® 7^7^)\E' = {"^lih ® 7^)^^, "^lih ® 
07^)^') = (pa, Ja), cx = 7*^7, with the wave functions of (6) and ^'o of (20) 
and compute pA- Inasmuch as in the chosen representation 7^ = — i7°7^7^7^ = 

- ° ) then we shall have pa = E •=i(i^j2 - iP,i)($°Vi$,). Under 

this situation for the ground state of P-meson described by (9) and (11) it is 
natural to put $j — ai^j which entails $°Vi<I>j — l/(47r) for any spinor of 
(11). Then we can define the weak axial form factor fA{K) for P-meson (by 
analogy with electric form factor in Section 3) as the Fourier transform of pa 



fA{K) = / pAe-'^'^d^x = E Q / r'^^'^-^^' 
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so computation along the same lines as for electric form factor of (15) yields 
the result 

f (K\=.A V (2/3,)"^-+V^r(a,) - Q,. sinKarctan(X//3,)] ^ 

^ (2/j,)"^+V^r(a,) (1 a] + 3a, + 2 > 

V3r(2a. + i) 'T^T^'WI^U ^""^ ^ 

(22) 

with Vj = Pj{l - gbj/Pj), Qj = Oj{l + gbj/Pj). It is clear from (22) that 
fA{K) is a function of K"^ and we can rewrite (22) for arbitrary 4- vector Q as 



= y (2/3,)"^-+V^r(a,-) - Q,- sinK-arctan(^|Q2|//3,)] 

(23) 

which passes on to (22) in the Breit frame where = —K^. 



5 Estimates for pcirameters of SU(3)-gluonic field in P- mesons 



5.1 Basic equations 



Now we are able to estimate parameters aj,bj,Bj of the confining SU(3)- 
field (3) in P-mesons. Though we can consider other possible functions of 
connected with the mentioned above current J a (on the analogy of the 
meson magnetic form factor, see Refs. [7,8]) but, obviously, the most natural 
function is /a(<5^) of (23). It is reasonable, therefore, to take for the function 
T{Q'^) mentioned in Subsection 4.1 just the function |/a(<5^)|, i-e., to put 
jr(Q2j ^ |/^(Q2)| so that we should have |/a(Q2 = ml)\ ^ 0.9364479961 
or 0.3236934271 respectively for P = tt^ and P = K^. Then, denoting the 
quantities mp/Pj = Xj, we obtain the following equation for parameters of the 
confining SU(3)-gluonic field in P-mesons 



\fA{Q'^ml)\^A, 



^ 2"i+^/^T{aj) Vj-Qj sin(ajarctana; 



0.9364479961, 
0.3236934271. 



(24) 



Finally, equations (12) and (16) should also be added to (24) and the system 
obtained in such a way should be solved compatibly if taking < r > ^ 0.672 
fm for P = TT^ and < r > ?a 0.560 fm for P = [9]. While computing for 
distinctness we take all eigenvalues Xj (j = 1, 2, 3) of the Euclidean Dirac 
operator T>o on the unit two-sphere equal to 1. 



13 



5.2 Numerical results 



The results of numerical compatible solving of equations (12), (16), and (24) 
are adduced in Tables 1-6 (with ji — mp). 

Table 1 

Gauge coupling constant, reduced mass jjlq and parameters of the confining SU(3)- 
gluonic field for 7r='=-mesons 



Particle 


9 


/xo (MeV) 


ai 


02 


6i (GeV) 


62 (GeV) 


Bi 


B2 


7r± ud. ikI 


().rM)l:-!l 


1.55172 


1), 017:501)2 


0, Oils 197 


0.178915 


-0.119290 


-0.2:-!0 


0.230 



Table 2 

Theoretical and experimental 7r^-meson mass and radius 



Particle 


Theoret. IJ, (MeV) 


Experim. ^ (MeV) 


Theoret. < r > (fm) 


Kxperim. < r > (fm) 


TT^ — ud, ud 


tj, = mu+mci + LOj (0, 0, 1) = 139.570 


139.56995 


0.673837 


0.672 



Table 3 

Theoretical and experimental values of the charged pion weak axial form factor and 
constant Aq 



Particle 


Theoret. |/a(Q^ = M^) 


Experim. |/a(Q^ = 


Aq (MeV^/^) 


TT^ — ud, ud 


0.9364898401 


0.9364479961 


12.72 



Table 4 

Gauge coupling constant, reduced mass jjlq and parameters of the confining SU(3)- 
gluonic field for ii'='=-mesons 



Particle 


3 


/xo (MeV) 


ai 


02 


61 (GeV) 


62 (GeV) 


Bi 


B2 


MS, us 


5.:-!()121 


2.20:-!87 


0.1()7182 


-0.0557501 


0.120150 


0.1:-!101(i 


-0.900 


0.290 



Table 5 

Theoretical and experimental i(r='=-meson mass and radius 



Particle 


Theoret. jj, (MeV) 


Experim. (MeV) 


Theoret. < T- > (fm) 


Experim. < r > (fm) 


US, US 


IX = ruu + rus + uJj{Q,0,l) = 493.677 


493.677 


0.544342 


0.560 



Table 6 

Theoretical and experimental values of the charged kaon weak axial form factor and 
constant Aq 



Particle 


Theoret. |/a(Q^ = M^)I 


Experim. |/a(Q^ = M^)| 


Ao (MeV^/2) 


— us, us 


0.3232652803 


0.3236934271 


2.88 
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6 Estimates of gluon concentrations, electric and magnetic colour 
field strengths 



Now let us remind that, according to Refs. [3,6], one can confront the field (3) 
with Too-component (volumetric energy density of the SU(3)-gluonic field) of 
the energy-momentum tensor (2) so that 

E'^ + H'' _ 1 f al + aia2 + aj 6^ + M2 + &i \ ^ A B 
"~ 2 ~2V ^ r^sinH J-r^^r^sinH 

(25) 

with electric E and magnetic H colour field strengths and with real A > 0, B > 
0. One can also introduce magnetic colour induction B — (Att x 10~^H/m) H, 
where H in A/m. 

To estimate the gluon concentrations we can employ (25) and, taking the 
quantity uj = F, the full decay width of a meson, for the characteristic fre- 
quency of gluons we obtain the sought characteristic concentration n in the 
form 

(26) 



n 



r 



so we can rewrite (25) in the form Too — T^^^ + Tq™ conforming to the con- 
tributions from the Coulomb and linear parts of the solution (3). This entails 
the corresponding split of n from (26) as n = ricoui + ^lin- 

The parameters of Tables 1 and 4 were employed when computing and for 
simplicity we put sin i? = 1 in (25). Also there were used the following present- 
day full decay widths - for 7r^-mesons: T = 1/r with the life time r = 2.6033 x 
10^*^ s and for A'^'^-mesons: F = 1/r with the life time r = 1.2386 x 10~^ s, 
whereas the Bohr radius ao = 0.529177249 • 10^ fm [9]. 

Tables 7-8 contain the numerical results for ricoui, JT-im, n, E, H, B for the 
mesons under discussion. 

Table 7 

Gluon concentrations, electric and magnetic colour field strengths in Tr^-mesons 



7r±— ud, ud: ro =< r >= 0.673837 fm 



r (fm) 



"coul (m ) 



"lin (m ) 



E (V/m) 



H (A/m) 



B (T) 



O.lro 



0.366804 X 10^^ 



0.141131 X 10^^ 0.380917 X 10^^ 0.201234 X 10^* 0.630981 X 10^^ 0.667250 X 10^^ 



0.366804 X 10^^ 



0.141131 X 10^^ 0.177811 X 10^^ 0.201234 X 10^^ 0.630981 X 10^° 0.667250 X 10^* 



0.756229 X 10° 



0.640814 X 10^^ 0.716437 X 10®^ 0.913716 X 10^^ 0.367794 X 10^° 0.449618 X 10^* 



lOr-o 



0.366804 X 10^* 



0.141131 X 10^1 0.141498 X 10^^ 0.201234 X 10^° 0.630981 X 10^^ 0.667250 X lO'^^ 



0.964380 X 10*^ 



0.228839 X 10^^ 0.228839 X 10^^ 0.326294 x 10^^ 0.676133 X 10^^ 0.849654 x lO'' 
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Table 8 



Gluon concentrations, electric and magnetic colour field strengths in i^^-mesons 



— us, us: 


ro =< r >= 0.544342 fm 




























r (fm) 


"coul 


"lin (m 


") 










E (V/m) 




H (A/m) 




B (T) 






O.lro 


0.303178 X 10^'' 


0.195700 


X 


10^^ 


0.305135 


X 


1067 


0.838745 X 


10=* 


0.906483 X 


10=1 


0.113912 


X 


10^6 


'■o 


0.303178 X 10^^ 


0.195700 


X 




0.498878 


X 


1063 


0.838745 X 


j(,22 


0.906483 X 


1020 


0.113912 


X 


10^5 


1.0 


0.266186 X 10^^ 


0.579876 


X 


j(,62 


0.846062 


X 


1062 


0.248527 X 


J022 


0.493437 X 


1(,20 


0.620071 


X 


10" 


lOro 


0.303178 X 10^^ 


0.195700 


X 


io«i 


0.198732 


X 


io«i 


0.838745 X 


1020 


0.906483 X 


lois 


0.113912 


X 


10" 




0.339454 X 10*^ 


0.207077 


X 


jq53 


0.207077 


X 


^q53 


0.887506 X 


ioi=^ 


0.932460 X 


IQl'i 


0.117176 


X 


10^° 



7 Chiral symmetry breeiking 



7. 1 Preliminaries 

As is known, in the late sixties of XX century, in light meson physics there arose 
notion of chiral symmetry and so far the latter has been actively exploited in 
phenomenology (see, e.g. reviews [15]). In its turn, to provide real world with 
chiral symmetry breaking there was supposed that chiral symmetry had been 
spontaneously broken. On the other hand, historically in fact simultaneously 
with notion of chiral symmetry there was created standard model (SM) of 
electroweak interactions which later quarks were also included in. We should 
note that SM (with one Higgs doublet) contains some description of chiral 
symmetry breaking: current masses of quarks (as well as lepton masses) are 
acquired through Higgs mechanism so for quark masses we obtain (without 
taking mixings into account) rriq = fgv/\^, where vacuum Higgs condensate 
V ~ 246 GeV. But little is known about coupling constants fq and much may 
be elucidated only with discovering Higgs bosons. 

As a result, we at present have two postulated mechanisms for spontaneous 
breaking of (global) chiral symmetry: one is associated with SM and another 
one should be related with the so-called Goldstone bosons accompanying vi- 
olation of any global symmetry. It is clear that both mechanisms should be 
connected in one or another way. As far as is known to us, up to now there 
has been not generally accepted (if any) recipe for reconciliation between men- 
tioned mechanisms. 

After QCD was created in the early seventies of XX century there appeared 
a hope that chiral symmetry breaking should be explained within framework 
of QCD and it should closely be related to confinement mechanism (see, e.g., 
Refs. [16]). Under the situation we can remark that our approach has a well- 
defined chiral limit with perfectly clear physical meaning and in view of this 
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it contains some description of chiral symmetry breaking. Let us look into it 
in more detail. 



1.2 Pion and kaon masses and chiral limit 



Let us more in detail write out expressions for pion and kaon masses from (12) 



^ , g'a.h , v/(l - gB,y - g'aj 

m^± =my, + md + — + iJ,o- rz —. 

1-gBi ll-^-Bil 



g'a2b2 , ^(l-gB^r-g^al 



g\ar + a2)(6i + 62) V + ^(^1 + ^2)]' " Q^^^i + 0.2? 

l + g{Bi + B2) \l + g{Bi + B2)\ 

murrid , . 

1^0 = , 27 

mu + md 



^ , g'a.b, , ^{l-gB,r-g^al 
mK± = m„ + ms + — + Hq- 



l-gB^ \l-gB^ 



mu + ms + — + Ho 



g^a2b2 J{l-gB2Y-g^ai 



l-gB2 1 - gB2 



g\ai + a2){bi + b2) ^J[1 + g{Bi + B2)f - g^ia^ + a2) 



mu + ms-\ — — „ . h Ho 



l + g{B, + B2) \l + g{Bi + B2)\ 

ixo = " . 28 

ruu + Vfis 

In chiral limit m„, m^, ms — > we obtain 

[m^^Urai ~ ~ ~ 1 + ^(Si + B2) ~ 130.8 MeV ^ , 

(29) 

/ X ^^«2&2 ^^(ai + a2)(6i + 62) oooA/rAA^n 

(m;..)...„. - - - 1 + ,(E, + E2) " ''^^ ^ ° 

(30) 

using the parameters g, aj, bj, Bj adduced in Tables 1 and 4 respectively. 

We can see that in chiral limit the pion and kaon masses are completely 
determined only by the parameters aj,bj,Bj of SU(3)-gluonic field between 
quarks, i.e. by interaction between quarks, and those masses have purely glu- 
onic nature ! As a consequence, if neglecting gluon field we exactly obtain 
{nip)chirai = 0. Analogously, the root-mean-square radii of pions (kaons) of 
(16) are well-defined in chiral limit and < r >chirai^ 0.673069 fm or 0.543223 
fm with parameters Oj, bj, Bj of Tables 1 and 4 respectively and those values 
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only slightly differ from 0.673837 fm or 0.544342 fm of Tables 2 and 5 accord- 
ingly at m„ , rriii , ^ 0. The same holds true for decay constants of (19) for 
leptonic decays fp = mp\fA{Q'^ = irip)] with = iTip)] of (24). So, even 

in chirally symmetric world the pions and kaons would have nonzero masses, 
the root-mean-square radii and decay constants fp for leptonic decays and all 
of those quantities would be determined only by SU(3)-gluonic interaction be- 
tween massless quarks, i.e. they would have a purely gluonic nature. It should 
be emphasized that we can neglect neither Coulomb electric colour field of so- 
lution (3) (parameters ai, 02) nor magnetic colour field (parameters 61, 62) or 
else effect does vanish according to (29)-(30), i.e. both parts of SU(3)-gluonic 
field of (3) are important for confinement and mass generation in chiral limit. 
Moreover, since gluons are verily relativistic particles then the most part of 
mass for light mesons is conditioned by relativistic effects, as is seen from (29) 
and (30). 

7.3 Quark condensate and chiral perturbation expansion 

Thus, our approach neatly says that chiral symmetry is rather rough approxi- 
mation holding true more or less only when neglecting SU(3)-gluonic field be- 
tween quarks. But we historically know that notion of chiral symmetry arose 
when no indications to existence of gluons were. So this was a reasonable ap- 
proach for its time. However, after discovering gluons in the late seventies of 
XX century, main constructions developed within chiral symmetry approach 
have not been essentially changed to take into account the fact of existence of 
gluons and such a situation is continuing up to now (see, e.g. reviews [15]). 

Let us consider, for example, how one should interpret one of the key relations 
of chiral symmetry approach - the Gell-Mann-Oakes-Renner relation connect- 
ing the so-called quark condensate < 0\uu + dd\0 >~ — 2(240 MeV)^ (an order 
parameter characterizing chiral symmetry breaking) with /^r, m„, rrid, (for 
more details see Refs. [15]) 



Then, under usual interpretation, since f^^O (experimental fact !) in chiral 
limit one should conclude that ~ and identify pion with a Goldstone 
boson. It is clear that absence of gluons was tacitly supposed because when 
deriving (31) no gluons were known about. 

As has been said above, however, presence of gluons between quarks entails 
that in chiral limit the left-hand side of (31) is well-defined and not equal to 
0. Under the circumstances 




'TT 



2 



{ruu + md) < Q\uu + dd\Q > . 



(31) 



< Oluu + ddlO >= 



(,fn^n^ chiral 

rUu + nid 



—00 , 



(32) 
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i.e., < 0\uu + dd\0 > becomes unphysical, unobserved parameter. Similar 
analogy: when the speed of hght c ^ oo relativistic mechanics passes on 
to Newtonian one where c is unobservable and no relations of Newtonian 
mechanics depend on c. Accordingly, the quark condensate is rather crude 
effective parameter that exists only when neglecting the fact of existence for 
gluons. We can, of course, try to amend (31) by changing + — >• m„ + 
md + {mT^)chirai with {mT^)chirai of (29) but all the same m.,^ will be ~ only if 
parameters of SU(3)-gluonic field a^, hj are very small. 

Also we can make a comment on the so-called chiral perturbation theory ( for 

more details sec reviews [15]). Within the latter approach, for example, mass 
square of P-meson (P = tt^, K"^) is sought in the form of chiral perturbation 
expansion 

m^P = xAp + yBp + 0{x'^,y'^) (33) 

with current masses x, y of quarks composing P-meson while coefficients Ap,Bp 
are calculated in accordance with the rather hazy rules [15]. We can, however, 
see from (27)-(28) that actual form of nip is 

ml^(A + x + y + ^B] ^ f{x,y), (34) 
V x + y ) 

where A, B depend only on gluonic field. But then standard differential cal- 
culus says to us that function /(x, y) (though being continuous at x = = 0) 
is not differentiable aX. x = y = Q and, accordingly, f{x,y) does not possess 
the Taylor expansion of the form (33) at the point x — y = 0. As a result, the 
expansion (33) is incorrect. 

We should, however, be somewhat careful: from nowhere it follows that pa- 
rameters aj,bj,Bj of the confining SU(3)-gluonic field in chirally symmetric 
world for pions and kaons should be the same as ones at present - the latter 
were evaluated in accordance with (12) at nonzero m„, m^, m^. Under the 
situation we can at least describe two scenarios of chiral symmetry breaking. 

7.4 First scenario 

In chirally symmetric world (e.g., in early stages of universe evolution) massless 
quarks interchange with gluons which generate nonzero masses of hadrons of 
purely gluonic nature, as was discussed above. After spontaneous breaking 
of symmetry in SM quarks acquire current masses through Higgs mechanism 
which entails additional contributions to hadron masses but parameters of 
gluon field between quarks remain the same. I.e., we suppose massless quarks 
to emit gluons in the same proportion as massive ones. As a result, parameters 
of gluon field in hadrons at present are the same as in chirally symmetric world. 
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1.5 Second scenario 

After spontaneous breaking of symmetry in SM massive quarks emit gluons in 
other proportion than massless ones and parameters of gluon field in hadrons 

at present are different from those in chirally symmetric world. But for to eval- 
uate the latter we should know, for example, pion mass in chirally symmetric 
world that is not equal to zero due to gluons as was discussed above. 

7.6 Concluding remarks 

At any rate, in either scenarios no additional mechanism of spontaneous sym- 
metry breaking connected with Goldstone bosons is required. Another matter 
that massless quarks differ from each other only by their flavours and we come 
to the problem of origin for flavours. It is clear, however, the latter problem 
cannot be resolved within QCD which just takes flavours as given from outside. 
So problem of origin of flavours requires coming out from QCD-framework. 
By the way, one possible solution of this problem from cosmological positions 
was proposed by us a long time ago [17]. 

To summarize, our confinement mechanism gives a physically reasonable ap- 
proach to problem of chiral symmetry breaking without any additional mech- 
anism of spontaneous symmetry breaking connected with Goldstone bosons. 



8 A possible relation with a phenomenological string-like picture 
of confinement 

8. 1 The confining potential and string tension 

The results obtained in Sections 5-6 allow us to shed some light on one more 
problem which has been touched upon in Ref . [6] . As is known, during a long 
time up to now there exists the so-called string-like picture of quark confine- 
ment but only at qualitative phenomenological level (see, e. g., Ref. [18]). Up 
to now, however, it is unknown as such a picture might be warranted from 
the point of view of QCD. Let us in short outline as our results (based on and 
derived from QCD-Lagrangian directly) naturally lead to possible justification 
of the mentioned construction. Thereto we note that one can calculate energy 
£ of gluon condensate conforming to solution (3) in a volume V through re- 
lation £ = JyTooT^ sin Mr dM(p with Too of (25) but one should take into 
account that classical Too has a singularity along z-axis (■»? = 0, tt) and we 
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Fig. 1. Vertical projection of region with the gluon condensate energy between 
quarks 

Fig. 2. Horizontal projection of region with the gluon condensate energy between 
quarks 

have to introduce some angle -i^o so -(^o < ''^ < tt — -i^o- As well as in Ref. [2], 
we may consider t^o to be a parameter determining some cone i) = -^q so the 
quark emits gluons outside of the cone. Now if there are two quarks Q11Q2 
and each of them emits gluons outside of its own cone '& — ■i?i^2 (see Figs. 1, 
2) then we have soft gluons (as mentioned in Section 1) in regions I, II and 
between quarks. 

Accordingly, we shall have some region V with gluon condensate between 
quarks Qi, Q2 and its vertical projection is shown in Fig. 1. Another projection 
of V onto a plane perpendicular to the one of Fig. 1 is sketched out in Fig. 2. 

Then, as is clear from Fig. 1, for distance R between quarks we have R — 
Ri sin i?i + i?2 sin -§2 and gluonic energy between quarks will be equal to 



V{R)= j Toqt'^ sinMrdMip = J J J (■:^ + -^^)^^^^drd'&dip+ 

V ri ^1 -ipi 

with constants A, B defined in (25). 

To clarify a physical meaning of the quantities ri^2 in Figs. 1, 2, let us recall 
an analogy with classical electrodynamics where is well known (see e. g. Rcf. 
[19]) that the notion of classical electromagnetic field (a photon condensate) 
generated by a charged particle is applicable only at distances much greater 
than the Compton wavelength Ac = 1/m for the given particle with mass m. 
Within the QCD framework the parameter Aqcd plays a similar part (see, 
e.g., Ref. [9,16]). Namely, the notion of classical SU(3)-gluonic field ( a gluon 
condensate) is not applicable at the distances much less than 1/Aqcd- In 
accordance with Section 3.1 we took Aqcd = A = 0.234 GeV which entails 
1/A ~ 0.8433 fm so one may consider ri^2 ^ O.lro where Tq is adduced in 
Tables 7-8 for charged pions and kaons. 

Under the circumstances, performing a simple integration in (35) with em- 
ploying the relations J di)/ sint? = lntan'i?/2, tan'i?/2 = sin'i?/(l + cost?) ~ 
(1 — cos'&)/ sin'j?, we shall without going into details (see also Ref. [2]) obtain 
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D^ 4:ipiA cos 1^ i l + cos^^j 
V(i?i,i?2) = Vo - 2^ ^ h 2^ 2^iBi?i In — -, (36) 

j^—l J^i j^—i -L COS i/j 

where Vo = E-=i Vo, = Eii (^^^-^ - 2(^,fir, In 'j^^). 

For the sake of simphcity let us put it!i = i?2, = '»?2 = '»?o, '/'i = </'2 = </?o- 
Then ^ ^ i?/(2sini?o) and from (36) it follows 

V{R) =Vo + ^ + kR (37) 
With a = -8(^o^sin2^o, = 2(^0-^; In 

We recognize the modeling confining potential in (37) which is often used when 
applying to meson and heavy quarkonia physics (see, e.g., Refs. [20,21]). Wc 
can, however, see that phenomenological parameters a, k, Vq of potential (37) 
are expressed through more fundamental parameters aj, bj connected with 
the unique exact solution (3) of Yang-Mills equations describing confinement. 
One can notice that the quantity k (string tension) is usually related to the 
so-called Rcggc slope a' = 1/ (27r/c) and in many if not all of the papers using 
potential approach it is accepted A; 0.18 GeV^ (see, e. g., Refs. [20,21]). 



8.2 Estimates of -Oq, (po for charged pions and kaons 



Under the situation we have the equation 

k = 2(^0-^ In 1±^^ « 0.18 GeV^ (38) 

Sm Vq 1 — COS Vo 

with B = {bl + 6162 + &2)/2, so let us employ (38) to estimate t^o; ^0 if using 
the results obtained in Table 1 and 4 for charged pions and kaons and also 
for the ground state of toponium rjt for that we use the parametrization from 
Ref. [5] with the values ai = 0.361253, as = 0.339442, bi = 48.9402 GeV, 
62 = 76.7974 GeV for the parameters of solution (3). Results of computations 
are presented in Table 9. 

If taking into account that only the values of ??0! fo between and 90° are 
of physical meaning and, according to Figs. 1, 2, the corresponding region 
V between quarks will be similar to a string-like one under the condition 
7''/2, <Po 0, then we can see from Table 9 that the characteristic 
transverse sizes D12 of the gluon condensate between quarks in fact tend to 
zero only in the case of heavy quarks, i.e., only for heavy quarks the gluon 
configuration between them might practically transform into a string. As a 
result, there arises the string- like picture of quark confinement but the latter 
seems to be warranted enough only for heavy quarks. It should be emphasized 
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Table 9 

Angular parameters determining the gluon condensate between quarks for charged 
pions, kaons and toponium ground state. 



Particle 


^0 




vr^ — ud, ud 


10° 


14.76° 




30° 


78.63° 




45° 


166.16° 




60° 


326.53° 


— us, us 


10° 


7.76° 




30° 


41.34° 




45° 


87.36° 




60° 


171.68° 


Vt — tt 


10° 


(0.305 X 10-4)° 




30° 


(0.162 X 10"^)° 




45° 


(0.343 X 10-3)° 




60° 


(0.675 X 10-3)° 




80° 


(0.240 X 10-2)° 




88° 


(0.123 X 10~^)° 



that string tension k of (38) is determined just by parameters 6i_2 of linear 
magnetic colour field from solution (3) which indirectly confirms the dominant 
role of the mentioned field for confinement. 

We cannot, however, speak about potential V{R) of (37) as describing some 
gluon configuration between quarks. It would be possible if the mentioned 

potential were a solution of Yang-Mills equations directly derived from QCD- 
Lagrangian since, from the QCD-point of view, any gluonic field should be 
a solution of Yang- Mills equations (as well as any electromagnetic field is by 
definition always a solution of Maxwell equations) . 

In reality, as was shown in Refs. [2,3] (see also Appendix C in Ref. [5]), po- 
tential of form (37) cannot be a solution of the Yang-Mills equations if si- 
multaneously a ^ 0,k 0. Therefore, it is impossible to obtain compatible 
solutions of the Yang-Mills-Dirac (Pauh, Schrodinger) system when inserting 
potential of form (37) into Dirac (Pauli, Schrodinger) equation. So, we draw 
the conclusion (mentioned as far back as in Refs. [4] and elaborated more 
in detail in Ref. [5]) that the potential approach seems to be inconsistent: 
it is not based on compatible nonperturbative solutions for the Dirac- Yang- 
Mills system derived from QCD-Lagrangian in contrast to our confinement 
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mechanism. Actually potential approach for heavy quarkonia has been his- 
torically modeled on positronium theory. In the latter case, however, one uses 
the unique modulo square integrable solutions of Dirac (Schrodinger) equation 
in the Coulomb field [condensate of huge number of (virtual) photons], i. e., 
one employs the unique compatible nonperturbative solutions of the Maxwell- 
Dirac (Schrodinger) system directly derived from QED-Lagrangian to describe 
positronium (or hydrogen atom) spectrum. 

To summarize, from the point of view of our approach both potential and 
string-like pictures of confinement arise only as some effective models derived 

in a certain way from the more fundamental theory based on exact solution 
(3) of SU (3)- Yang-Mills equations. This conclusion is in concordance with the 
preliminary one obtained in Ref. [6]. 



9 Discussion and concluding renicirks 

9. 1 Discussion 

As is seen from Tables 7-8, at the characteristic scales of charged pions and 
kaons the gluon concentrations are huge and the corresponding fields (electric 
and magnetic colour ones) can be considered to be the classical ones with 
enormous strengths. The part ncoui of gluon concentration n connected with 
the Coulomb electric colour field is decreasing faster than nun, the part of n 
related to the linear magnetic colour field, and at large distances niij, becomes 
dominant. It should be emphasized that in fact the gluon concentrations are 
much greater than the estimates given in Tables 7-8 because the latter are 
the estimates for maximal possible gluon frequencies, i.e. for maximal possible 
gluon impulses (under the concrete situation of charged pions and kaons) . As 
was mentioned in Section 1, the overwhelming majority of gluons between 
quarks is soft, i. e., with frequencies much less than F = 1/r ^ 0.253 x 10~^ 
eV, for example, in the case of pions, so the corresponding concentrations are 
much greater than those in Tables 7-8. The given picture is in concordance 
with the one obtained in Refs. [4-8]. As a result, the confinement mechanism 
developed in Refs. [1-3] and described early in Section 1 is also confirmed by 
the considerations of the present paper. 

It should be noted, however, that our results are of a preliminary character 
which is readily apparent, for example, from that the current quark masses (as 
well as the gauge coupling constant g) used in computation are known only 
within the certain limits and we can expect similar limits for the magnitudes 
discussed in the paper so it is necessary further specification of the parameters 
for the confining SU(3)-gluonic field in charged pions and kaons which can be 
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obtained, for instance, by calculating the width of decay tt^ 7r° + e^ + t'e(t'e) 
with the help of wave function of 7r*^-meson discussed in Ref. [7] and so on. 
We hope to continue analysing the given problems elsewhere. 

9.2 Concluding remarks 

The results of present paper as well as the ones of Refs. [4-8] allow one to speak 
about that the confinement mechanism elaborated in Refs. [1-3] gives new 
possibilities for considering many old problems of hadronic (meson) physics 
(such as nonperturbativc computation of decay constants, masses and radii 
of mesons, chiral symmetry breaking and so forth) from the first principles of 
QCD immediately appealing to the quark and gluonic degrees of freedom. This 
is possible because the given confinement mechanism is based on the unique 
family of compatible nonperturbativc solutions for the Dirac- Yang- Mills sys- 
tem directly derived from QCD-Lagrangian and, as a result, the approach 
is itself nonperturbativc, relativistic from the outset, admits self-consistent 
nonrelativistic limit and may be employed for any meson (quarkonium) . 

The given paper to a certain degree summarizes studying nonet of pseudoscalar 
mesons realized in Refs. [6-8] within the framework of our approach and we 
can ascertain the fact that, on the whole, this nonet can be described from the 
unified point of view of our confinement mechanism. In line with the above, 
obviously, one should now pass on to vector mesons (p, 0, uj...) and also to 
the light scalar mesons whose nature has been controversial over 30 years [22]. 
As is clear from Section 2 and Appendices A, B, there exists a large number 
of relativistic bound states in the confining SU(3)-gluonic field (3) so all the 
mentioned mesons can probably correspond to some of those states and be 
described by their own sets of parameters aj, bj, Bj of solution (3). Finally, 
one should think about possible ways to extend the approach over baryons, in 
particular, over nucleons. 



Appendix A 

Wc here represent some results about cigenspinors of the Euclidean Dirac 
operator on two-sphere S'^ employed in the main part of the paper. 

When separating variables in the Dirac equation (4) there naturally arises 
the Euclidean Dirac operator Vq on the unit two-dimensional sphere §^ and 
we should know its eigenvalues with the corresponding cigenspinors. Such a 
problem also arises in the black hole theory while describing the so-called 
twisted spinors on Schwarzschild and Reissner-Nordstrom black holes and it 
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was analysed in Rcfs. [3,23], so we can use the results obtained therein for our 
aims. Let us adduce the necessary relations. 

The eigenvalue equation for corresponding spinors $ may look as follows 

Vo^ = A$. (A.l) 

As was discussed in Refs. [23], the natural form of (arising within applica- 
tions) in local coordinates i?, (p on the unit sphere §^ looks as 



1028^ + i(T^^^ (d^ - ^(J2(73COS^9 

sm^y V 2 



„ 1 _ cot / A ^\ 

(7i(72d^ + - — ^CTicrsO^ (A.2) 

smv 2 



with the ordinary Pauli matrices 



on /o-A /lo 



so that aiVo — — "DqCi- 

The equation (A.l) was explored in Refs. [23]. Spectrum of Dq consists of the 
numbers A — ±(/-|-l) with multiplicity 2(/ + l) of each one, where 1 — 0,1, 2, .... 
Let us introduce the number m such that — / < m < /+! and the corresponding 
number m' = m — 1/2 so \m'\ < I + 1/2. Then the conforming eigenspinors of 
operator Vq are 

*=(*:)= '■t'^ = f ( 5'""' 1 5'"^ ) e""'- (A.3) 

U+1 



with the coefficient C = y and k — 1 + 1/2. These spinors form an or- 
thonormal basis in LKs^) and are subject to the normalization condition 

TT 27r 

J^^MQ^J y^(|$i|^ + |$2|^)sin«d(^= 1. (A.4) 

§2 

Further, owing to the relation aiVg — —T>q(T\ we, obviously, have 

cri^^\ = ■ (A.5) 



As to functions P^/„/(cos'i?) = P^, ^,{cos'd) then they can be chosen by mis- 
cellaneous ways, for instance, as follows (see, e. g., Ref. [24]) 



^m'n'(cOS^) 



{k - m')\{k - n')\ (l + cos-d^ 
\ {k + ■wl)\{k + niy. \ l-cos'&^ 



m'+n' 
2 



X 
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.^.ir-w)(^-i)'(i--')!(i--')!V 2 ; ^ • ^ 

with the orthogonahty relation at m', n' fixed 

TT 

/ P;t(cos^)P^',„,(cos^) sin« = ^^S,,, . (A.7) 



It should be noted that square of T>q is 

COS 1^ 11 
Vl = -As272 + (72^3-^9^ + — ^ + - , (A.8) 
sm XT 4sin 17 4 

while laplacian on the unit sphere is 



1 



sin + —^d' =di + cot + —^d' , (A.9) 



^ sin^? " sin'^ ^ " sin^i? 

so the relation (A.8) is a particular case of the so-called Weitzenbock-Lichnerowicz 
formulas (see Refs. [25]). Then from (A.l) it follows = and, when 

using the ansatz $ = P(T?)e-^"''<^ = (p ) e"'"*''^, Pi,2 = PiA'^), the equation 

= turns into 

/ o2 no + i m' cos -(9 \ „ 

V sm^i? sm^^? y 

A^-i)p, (A.IO) 
wherefrom all the above results concerning spectrum of Vq can be derived [23] . 

When calculating the functions P^/„/(cos'i?) directly, to our mind, it is the 
most convenient to use the integral expression [24] 



^m'n'(cOS^) = ^ 



{k - m')\{k + m')] /•/,,2 ^ . -W2 . A 
27r^ {k - n')\{k + n')\ J \ 2 2 J 



and the symmetry relations {z — cos •&) 



k+n' 



im ip 



dip 



(A.ll) 



Pm'n'i^) — Pn'm'i^)^ Pm'-n'i^) — P-m'n'i^)^ Pm'n'i^) — P-m'-n'i^) > 



pk ( „\ _ •2fe-2m'-2n' pfc / \ 



(A.12) 
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In particular 



2m 



{2k)\ 



— sin"-'' {'d/2) cos^'+"' {'d/2) . (A. 13) 



\ {k - n')\{k + n')\ 



Eigenspinors with A = ±1, ±2 



If A = ±(/ + 1) = ±1 then Z = and from (A.3) it follows that k^ 1+1/2^ 

1 /2 

1/2, \m'\ < 1/2 and we need the functions P^, _j_^y2 ^^^^ easily evaluated 
with the help of (A.11)-(A.13) so the eigenspinors for A = — 1 are 



^ ^ C / cos ^ + i sin ^ X ^ ^ ^ / cos - + ^ sin - \ ^_,,/, 

2 V cos I - z sm I y ' 2 V - cos I + ? sm I y ' ^ ^ 

while for A = 1 the conforming spinors are 



^ ^ ^ cos I - i sin I N ^,,/, ^ ^ C; / - cos f + . sin | j 

2 V cos I + z sm I y ' 2 V cos I + z sm I ' ^ ^ 



with the coefficient C = ^Ji/{2t:). 

It is clear that (A.14)-(A.15) can be rewritten in the form 

A = -1 : ^ = ? f 't. 1 e^^'\ or $ = ^ f \ e-^-'\ 



so the relation (A. 5) is easily verified at A = ±1. 

In studying vector mesons and excited states of heavy quarkonia eigenspinors 
with A = ±2 may also be useful. Then k = I + 1/2 = 3/2, \m'\ < 3/2 and 
we need the functions P'^^'j^^i^ ^^^^ ^an be evaluated with the help of (A.ll)- 
(A.13). Computation gives rise to 

3/2 . Q . 'i? 3/2 

-^3/2,-1/2 ^ -— sin'j^sin - = -P_3/2,i/2' 

p3/2 _ V3 • n _ rfill 

-^3/2,1/2 ~ 2 ^^^''^COS 2 — -r_3/2_i/2' 

d3/2 % ( . ■& o • 3 \ 3/2 

^1/2,-1/2 ^ ~4 2 ~ 2 J " P -1/2,1/2^ 
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<2,i/2 = \ (cos ^ + 3 COS = P!f/2,-i/2, (A.17) 
and according to (A. 3) this entails eigenspinors with A = 2 in the form 

- f -^S^ - e-^-/^ sin ^ f ] e-^^^/^ (A.18) 

8 V 3e*~ + e-*2 y '22 \ e'2 J ^ ^ 

with C = l/^/vT; while eigenspinors with A = — 2 are obtained in accordance 
with relation (A. 5). 



Appendix B 



We here adduce the explicit form for the radial parts of meson wave functions 
from (6). At rij — they arc given by 



while at rij > by 



|j,F,2 = zC,g,r"^e-^^- (1 + 1 



(B.l) 



Fj2 = iCjQjr'^^e 



(B.2) 



with the Laguerre polynomials LP^(rj), rj — 2(3 (3j — jiQ — u'j -\- at 

j = 1, 2, 3 with 63 = ~{bi + 62), Pj = gbj + Pj, Qj = A*o - ^j, Yj = PjQj<^j + 
(P/ - Qj^)^aj/2, = P.gjAj- + (P/ + Qj)gaj/2 with 03 = -(ai + as), A^ = 

\j - gBj with P3 = -(Pi + B2), oij = ^A| - g'^aj, while A^- = ±{lj + 1) 
are the eigenvalues of Euclidean Dirac operator on unit two-sphere with 
Ij = 0,1,2,... (see Appendix A) and quantum numbers rij = 0,1,2,... are 
defined by the relations 

^ ^ gbjZj - (3jYj ^ ^^^^ 
PjPjQj 

which entails the quadratic equation (7) and spectrum (8). Further, Cj of 
(B.1)-(B.2) should be determined from the normalization condition 

00 

li\Fji\' + \F,2\')dr^-. (B.4) 
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As a consequence, we shall gain that in (6) G L'^M.''') at any t G M and, 
accordingly, \E' = (^1,^2! ^3) niay describe relativistic bound states in the 
field (3) with the energy spectrum (8). As is clear from (B.3) at rij = we 
have ghjl Pj = Yj/Zj so the radial parts of (B.l) can be rewritten as 

More details can be found in Refs. [1,3]. 
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